Electron self-trapping and fluctuation density-of-states tail at the critical point 
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We consider electron self-trapping due to its interaction with order-parameter fluctuations at the 
second-order phase-transition or critical point (for example, at the Curie temperature in magnetic 
or ferroelectric semiconductors) . Using Feynman path integral approach the autolocalization energy 
and the size of the self-trapped state (fluctuon) are estimated. It is shown that the fluctuon states 
are connected with the Lifshitz tail of the electron density-of-states, the parameters of this tail being 
determined by the critical exponents. 
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The interaction of the charge carrier in a semiconduc- 
tor with some order-parameter fluctuations can drasti- 
cally change its state leading to a self-trapping, or au- 
tolocalization 0, S 0> IS S • This phenomenon is 
of crucial importance, for example, for so hot subject as 
phase separation in magnetic semiconductors and colos- 
sal magnetoresistance materials 0, 0- , where the mag- 
netization plays the role of the order parameter. Since 
the band motion of the electron is easier (and, hence, the 
bandwidth is larger) for ferromagnetically ordered state 
the electron in antiferromagnetic or magnetically disor- 
dered surrounding creates a ferromagnetic region (mag- 
netic polaron, "ferron" 6j, or "fluctuon" 0]) and turns 
our to be self-trapped in this region. Recently, the for- 
mation of the magnetic polaron in ferromagnetic semi- 
conductors EuB 6 and, possibly, EuO, has been observed 
(see Refs0,0 and references therein). The order param- 
eter can be also of different origin, e.g., electric polariza- 
tion in ferroelectric semiconductors, or crystallographic 
order parameter in ordered alloys |2|. Phase transitions 
in helium monolayers |Toj and ultracold atom systems in 
optical lattices 11] might be novel interesting examples. 

The "driving force" of the self-trapping is always a 
band narrowing in disordered state in comparison with 
the completely ordered one. This intuitive view can 
be confirmed by a rigorous consideration for the classi- 
cal s — d exchange ( "Kondo lattice" ) model on the Bethe 
lattice however, for the real lattices it is actually not 
very accurate. It was demonstrated already in a semi- 
nal paper by Brinkman and Rice Q that the band edges 
do not depend on the degree of spin disorder and the 
band narrowing means rather the decrease of the elec- 
tron density-of-states (DOS) momenta, the DOS near the 
edges being exponentially small in the disordered state 
(Lifshitz tails 0, Q] ) . The fluctuons (we will use this 
term to emphasize that the magnetic origin of the order 
parameter is actually not relevant for the problem under 
consideration) are connected with the Lifshitz tails. 

The autolocalization region becomes larger when the 
energy goes closer to the band edge (and the edge it- 
self corresponds to the state with complete ordering in 



the whole crystal) Q. The consideration of the fluc- 
tuon problem is crucially dependent on the ratio of the 
autolocalization radius I to another scale of the system, 
namely, the correlation length £. In particular, a stan- 
dard phenomenological consideration 2j is inapplicable 
at the critical point (or the point of the second-order 
phase transition), where £ — > oo. At the same time, this 
case is especially interesting since the fluctuations near 
the critical point are strongest and the self-trapping con- 
ditions fulfill much easier in this case. Here we present a 
theory for the Lifshitz DOS tails and fluctuon states at 
the critical point T = T C . 

Following Refs. 0,0 we will use the path integral varia- 
tional approach developed by Feynman for the problem of 
polaron in ionic crystal [bsl Il6| . For simplicity, we start 
with the case of a scalar order-parameter acting only on 
the orbital motion of the electron and not on its spin (for 
example it may be the critical point in ferroelectrics); 
some generalization will be considered further. Then, in 
continuum approximation, the Hamiltonian of the sys- 
tem consisting of the electron and the order-parameter 
field can be written in a simple form 

H f (ip)+H e {r,<p), H e (r,tp) = ' 
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9<P(*) (1) 
. m is the 



where we have chosen the units h = m 
electron effective mass, r is the electron coordinate, tp (r) 
is the order-parameter field with its own Hamiltonian 
Hf (ip) and g is the coupling constant; we will choose it in 
such a way that tp varies between —1 and 1. The partition 
function of the whole system may be transformed to 
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where (3 — T 1 is the inverse temperature, Zf = 
Tr (/ ,e~ /3W ^ ¥ '' is the partition function of the field, 
tp (r,r) = z TH tW>tp (r) e - TH f^ and 



(A(<p)) f = ^Tr v e-^f^A(<p) 
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is the average over the field states. Here we will consider 
only classical case where r-dependence of the field can be 
ne glec ted. Using Feynman path-integral approach 0, 
Il7l Il8| and taking average over (p yields for the electron- 
only free energy 

T = -~{\nZ-\nZf) = -~\n J e~ s V [r (r)] , 



r(0)=r(/3) 

where Sq + Si„t is the effective action, 
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IC m (r (n) ...r (r m )) d,T 1 ...d,T„ 



and /C m (ri; r m ) are the m-th cumulant correlators, 
defined recursively by 

/Ci (n) = (r!)) /; 
IC 2 (ri;r 2 ) = (n) p (r 2 )) / - /Ci (r x ) K x (r 2 ) , ... (6) 

etc. Further we consider only the case of JC\ = 0. 

To estimate J- we use the same trial action as in Refs0, 
4, S t — S + Spot where 

S P ot = ^J^ j\r{T)-r{o)] 2 drda, (7) 

the oscillator freq uency u> being trial parameter (in con- 
trast with Refll5l we do not introduce any retardation in 
the trial action since our field is supposed to be static). 
Then the Peierls-Feynman-Bogoliubov inequality applied 
to Eq. (QJ reads 



F < J~ t + — (Si n t — S po t) t 
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where 



F t = -\\n( e -^P[r(r)], 

P Jr(O)=r(0) 
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A[r (T)]e^- 5f P[r (r)] 

r(0)=r(/3) 

which is equivalent to 
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To proceed, we will pass to the Fourier transforms of 
the cumulants IC m (Ki, ..,K m _i) and take into account 



that for the Gaussian trial action St one has 
^exp J i ^2 K j I r ( t j) ~ r ( T ™)] 



exp 
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where D is the space dimensionality, are the wave- 
vectors. 

For the states in the tail the variational parameter u> 
satisfies the inequalities (3uj 3> 1 and u> « W — 5^,n ax 
where K max is the Debye wave vector and W is of order 
of the electron bandwidth (the last inequality just means 
that the fluctuon size I = l/y2uj is much larger than 
interatomic distance) . Thus we have @, Q 
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where is the unit lattice volume. If one can ne- 
glect the K-dependence of the cumulants this sum can 
be transformed 0, to the phenomenological expression 
in terms of the fluctuation free energy in homogeneous 
field created by the electron. This assumption works not 
too close to the critical point where I >■ £. At the crit- 
ical point this assumption can be never used. Instead, 
the scaling properties hold 

K m (Ki, ...,K m _!) = a( 2 -")( m - 1 )/C m (oKi, aK m ^) 

(13) 

where 77 is the anomalous dimensionality critical expo- 
nent. Making replacement of variables by K 3 = Jq>Cj, 
where q = lu/W, and further returning to the real space 
coordinates conjugate to ttj, we get 
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where d = D — 2 + 77 is the anomalous space dimension- 
ality, 
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is the potential localized within Debye sphere around 
r = 0; f[U(r)] is the change of the free energy of the 
fluctuations bath upon switching an external potential 
U(r) on. 

Earlier |j| we considered the Gaussian approximation 
(only the term with m = 2 in Ea.l|14|)'l which is applicable 
for not too large coupling constants. To consider the 
states near the band edge E — —g we have to sum up the 
series ifT^jl . To this end, we employ the fact that Un (r) is 
virtually D-dimensional delta function for typical r > I. 
Further, to avoid ultraviolet divergencies, we return to 
the discrete-lattice Ising model, where ip ri = ±1. In this 
framework we adopt that the above potential acts only 
at site Ti = 0, viz. ud (r,j) = g<5 rii o- This assumption 
can only affect some numerical coefficients of order of 
unity in the following estimations. Then in the limit of 
gqd/? i corresponding to the states near the band edge 
one obtains 



/ q d/2 u D (r) 
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where Gd = (<^(o),i)^ • This is just a number larger 
than 1 which can be calculated if all necessary correlation 
functions are known. For the Ising model with nearest- 
neighbor interaction J, one can easily derive 

Gn 1 = exp {(3 c T f ) (j[ (cosh/3 c J + tp ri sinh/3 c J)^ 

(17) 

where z is the nearest-neighbor number and J- / is the 
Ising model free energy per site at T = T c = f3~ . 

Substituting Eq. (| 1 fcil) into Eq. I|14fl and minimizing with 
respect to q one finds for the optimal fluctuon size 
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and for the fluctuon energy 
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The fluctuon forms only at Tq < 0. This requirement 
yields necessary condition for the self-trapping, which 
reads 
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where the band is assumed wide in the sense that W ^S> g. 
In the case opposite to one given by Eq. lt2"UI) the Gaussian 
approximation works 

The estimation ()16(l does not work for the systems with 
continuous broken symmetry such as XY or Heisenberg 
model; in that case some logarithmic corrections arise. 
To be specific, let us consider the case of the XY model 



where ip — cos 8 with the angle 9 distributed on the in- 
terval [0,27r). Then, instead of Ea. l|15|l . the following 
asymptotic takes place 
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(S (6»o )) f , where 



with a constant Gd given by G~jj — vvuj/t 
the average is with XY model on discrete lattice. After 
minimization one has, instead of Ea. i|19|) 
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This expression can be applied, for example, for the 
case of two-dimensional XY model up to the Kosterlitz- 
Thouless transition [l^. Eof where the correlation func- 
tions decay by power-low with the distance and d = r\ 
grows linearly with the temperature. 

Using the results obtained for the fluctuon energy 
one can restore the asymptotic of DOS N (E) near the 
band edge E — —g. To this end we use the Laplace 
transformation connecting the partition function Z (/?) = 
exp \-&T (/?)] and the DOS 



Z(f3) = / N(E)e 



dE, 
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inverse transformation, and the saddle point method; the 
corresponding asymptotics are connected by the so called 
Tauberian theorems; this approach was used by Fried- 
berg and Luttinger [2l| to obtain the Lifshitz tail for 
disordered systems. To employ this in our treatment, we 
consider the inverse temperature j3 in Eqs. (|19|) . 122|) 
as a parameter, except the constant Gd and, possibly, 
critical exponents. Acting similarly to the derivation of 
Eq.(2.19) of Ref.jiJ] one can find the Lifshitz tail in the 
critical point for the discrete order parameter 



N (E) oc exp 
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where in both cases e = E + g. 

Our final results mean that for the case of discrete 
order parameter such as in the Ising model the main dif- 
ference between the Lifshitz tails at the critical point and 
above the critical region is just a replacement of the space 
dimensionality D by the anomalous space dimensionality 
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d = D — 2 + r\. For the case of continuous order parame- 
ter additional logarithmic stretching of Lifshitz exponent 
emerge. 
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